Anomalous enhancement of tetragonality in PbTiOs induced by negative pressure 
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Using a first-principles approach based on density-functional theory, we find that a large tetragonal 
strain can be induced in PbTi03 by application of a negative hydrostatic pressure. The structural 
parameters and the dielectric and dynamical properties are found to change abruptly near a crossover 
pressure, displaying a "kinky" behavior suggestive of proximity to a phase transition. Analogous 
calculations for BaTiOs show that the same effect is also present there, but at much higher negative 
pressure. We investigate this unexpected behavior of PbTi03 and discuss an interpretation involving 
a phenomenological description in terms of a reduced set of relevant degrees of freedom. 
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I. INTRODUCTION 

Recent work has shown that single-crystal solid- 
solution ferroelectric perovskites can have dramati- 
cally improved electromechanical properties compared 
to conventional transducer materials Representa- 
tive materials include PbZn 1 / 3 Nb 2 /303-PbTi03 and 
PbMg!/ 3 Nb2/303-PbTi03, which have ultrahigh piezo- 
electric coefficients and low dielectric loss. These ma- 
terials have also been observed to exhibit electric-field 
induced phase transformations to "ultrahigh" strain 
statesi 

PbTi03 serves as the common parent compound for 
this class of materials, and may be supposed to play 
an important role in the observed behavior. Since the 
discovery of ferroelectricity in perovskite oxides in the 
1950's, PbTi03 has been the focus of extensive experi- 
mental and theoretical study, ft has a single phase tran- 
sition at T c = 766 K from a paraelectric cubic phase 
to a ferroelectric tetragonal phase, and a c/a of 1.06 
at low temperature. The structure and properties of 
PbTiC>3 have been widely studied using first-principles 
calculations AS&LSi Nevertheless, we report here a fea- 
ture of the behavior of PbTiC>3 that had not previously 
been noticed. Our calculations show that an enormous 
tetragonal strain can be induced in PbTiC>3 by applica- 
tion of a negative hydrostatic pressure. The structural 
parameters, such as cell volume and atomic displace- 
ments, are found to change abruptly near a crossover 
pressure, displaying a "kinky" behavior suggestive of 
proximity to a phase transition. Analogous calculations 
for BaTiC>3 show that the same effect is also present 
there, but at much higher negative pressure. 

In this paper, we investigate this unexpected behavior 
of PbTiC>3, and discuss its interpretation using a phe- 
nomenological description in terms of a reduced set of 
relevant degrees of freedom. We make the notion of prox- 
imity to a phase transition more precise by demonstrat- 
ing that small changes in the parameters in this descrip- 
tion can take the system through a triple point, leading 
to first-order transition behavior. Although the applica- 
tion of negative pressure is not feasible experimentally, 
our theoretical study provides useful insights into the 



structural instabilities of PbTiC>3, and may ultimately 
help suggest other, more practical avenues leading to en- 
hanced tetragonality in PbTiC>3 and related compounds. 

The paper is organized as follows. Section II provides 
the technical details of our first- principles calculations. 
Sec. Ill gives the results of our computations for the 
tetragonal distortion, structural parameters, and lattice 
dynamical properties as a function of applied negative 
pressure. In Sec. IV we introduce a phenomenological 
model that allows us to identify the relevant degrees of 
freedom responsible for this behavior, and we conclude 
with a brief discussion in Sec. V. 



II. COMPUTATIONAL DETAILS 



All ab initio calculations were performed in the frame- 
work of the Hohenberg-Kohn-Sham density-functional 
theory (DFT) within the local-density approxima- 
tion (LDA). We use the ABIIMIT package^ a plane- 
wave pseudopotential code that, in addition to 
ground-state total-energy and force calculations, al- 
lows linear-response computations of phonon frequen- 
cies and Born effective charges. Our calculations use 
the Perdew-Wang 10 parameterization of the Ceperley- 
Alder— exchange-correlation, and the extended norm- 
conserving pseudopotentials of Teter>i£ These pseudopo- 
tentials include the O 2s and 2p, the Ti 3s, 3p, 3d and 
4s, the Ba 5s, bp and 6s, and the Pb bd, 6s and 6p in 
the valence states. We have used an energy cutoff of 
60 Ha throughout. The integrals over the Brillouin zone 
have been replaced by a sum over a 6x6x6 k-point mesh. 
Convergence of the relaxations requires the Hellmann- 
Feynman forces to be less than 0.003 eV/A. We have 
computed the eigenvalues and eigenvectors of zone-center 
force-constant matrices by using both finite-difference 
(frozen-phonon) and linear-response approaches (with 
typical displacements of ±0.007 a. u. for the former), find- 
ing excellent agreement between the two schemes. 
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FIG. 1: Negative-pressure dependences of (a) volume (in 
Bohr 3 ); (b) tetragonal strain 773 and 771; (c) atomic displace- 
ments in the z direction (in c units) for tetragonal PbTiOs. 



III. RESULTS 
A. Structural response 

The observed crystal structure of ferroelectric PbTi03 
has space group PAmm, with Pb in the (la) Wyckoff 
position (0, 0, Ti in (lb)(|, |, | + £ 2 ), and O in 
(lb)(i,i,£ 3 ) and (2c) (i,0,i+£ 4 ) (0, ±,±+£4). The 
free structural parameters are the a lattice constant, 
the c/a ratio and the atomic displacements along z, £i 
(expressed in units of c). At zero pressure, our LDA 
calculation yields a T=Q equilibrium lattice constant of 
a = 7.301 a.u. and c/a — 1.037, which are ~ —1% and ~ 
—3% less than the experimental values of 7.373 a.u. and 
1.065, 13 > respectively. The fully-relaxed internal coordi- 
nates are £ Pb = 0.0579, f Ti = 0.0268, £_ 0l ,o 2 = -0.0335 
and £o 3 = —0.0177, where Oi and O2 are the "in-plane" 
oxygens (2c) and O3 is the apical (along 2, (lb)) oxygen of 
the Ti-centered oxygen octahedron. We use the conven- 
tion that ^2 k £fc = 0. Throughout this work the reference 
state will be our theoretical minimum-energy ideal cubic 
structure (ao = 7.331 a.u.), in terms of which the strains 
are defined as rji — (a — ao)/ao and 773 = (c — ao)/arj. 

We perform full, unconstrained optimization of the 
structural parameters of tetragonal PbTiC>3 as a func- 



tion of external pressure ranging from to —7 GPa. At 
a given pressure, the set of parameters that minimizes the 
total energy is determined. The negative-pressure depen- 
dence of the optimized structural parameters is shown in 
Fig. 2] where panels (a)-(c) display the unit-cell volume, 
strains (771 and 773), and internal atomic displacements in 
the z direction, respectively. As can be seen from the 
figure, all of the structural parameters display an abrupt 
change around a crossover pressure p c ~ —4.8 GPa in a 
way that suggests proximity to a phase transition. The 
rapid enlargement of the unit-cell volume results from 
an abrupt increase of 773, while it should be noted that 
the in-plane strain 771 decreases slightly in the same pres- 
sure range. The net effect of the negative pressure is to 
stretch the unit cell strongly along z and slightly squeeze 
it in the plane. For example, just above the transition 
at p = -5.3 GPa, we find 771 = -0.02 and 7/3 = 0.186, 
and the resulting c/a is 1.21. Next, considering the re- 
laxed atomic positions, the most remarkable feature is 
the change in character of the oxygen-displacement pat- 
tern. At pressures below p c the apical oxygen is dis- 
placed less than the in-plane oxygens, whereas at pres- 
sures above p c it is displaced more. At p c both types of 
oxygens are displaced by nearly the same amount, which 
means that the oxygen cage forms a tetragonally strained 
octahedron. 

For each of the optimized structures, we have also 
carried out spontaneous polarization calculations (up- 
per panel of Fig. [2J as a function of negative pressure 
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FIG. 3: Negative-pressure dependence of volume (in Bohr 3 ) 
for tetragonal BaTi0 3 and PbTi0 3 . 



using the Berry-phase theory of polarization^* It can 
be seen that the polarization shows the same kinky be- 
havior at p c . To get some insight into the polarization, 
we compute the Born effective charges (Z£). The bot- 
tom panel of Fig. [3 shows Z£ 33 in the relaxed tetrago- 
nal structure at different pressures, as obtained by us- 
ing density-functional perturbation theory. The results 
at zero pressure for k = Pb, Ti, Oi and 3 are 3.84, 
6.77, —2.55, —5.50, similar to previously calculated val- 
ues of 3.92, 6.71, —2.56, — 5.51 ji* respectively. However, 
a direct comparison is not possible because the latter val- 
ues were computed at the experimental tetragonal lattice 
constants. The most important features are the anoma- 
lously large effective charges of Ti and O3 (oxygen along 
the bond) compared with their nominal charges (+4 and 
—2 respectively) and the anisotropy of the oxygen charge. 
These features persist throughout the negative-pressure 
range of interest, although all of the effective charges ap- 
proach somewhat closer to the nominal values above p c . 
Most notably, Z^ i and Zq 3 change by ~ 25% while pass- 
ing through p Cl suggesting a weakening of the Ti-0 bond, 
whereas Z* Ph and Zq i decrease less noticeably. This 
follows the same trend observed in other perovskites, 16 
where the Z*'s decrease as the ions are displaced away 
from their high-symmetry cubic sites. The large changes 
in Z^ ri and Zq 3 reflect the change of the Ti environment 
along the Ti-0 chains; the shortened T1-O3 bonds re- 
main almost constant in length while passing through p c , 
while the elongated bonds lengthen abruptly there. Note 
that even though the Z*'s decrease, the enhancement of 
the ionic displacements is so strong that the overall effect 
is a marked enhancement of the spontaneous polarization 
with negative pressure. 

We have also investigated BaTiC>3 under negative pres- 
sure. In order to facilitate comparison with the re- 
sults for PbTiC>3, we have imposed tetragonal symmetry, 
even though the ground state of BaTiC>3 is rhombohe- 
dral. Interestingly, we find the same anomalous effect 



as in PbTiC>3, but it occurs at much higher pressure in 
BaTiC>3. The computed optimized volume of BaTiC>3 in 
the tetragonal phase is plotted as a function of negative 
pressure in Figure |3 together with the corresponding re- 
sults for PbTiC>3 for comparison. As seen, in BaTiC>3 the 
jump of the volume occurs at p c ~ —10.6 GPa, a pressure 
that is approximately twice as large as for PbTiOs. Fur- 
thermore, the abrupt volume enhancement is even larger; 
in a pressure interval of width 1 GPa around p Cl the vol- 
ume jump is ~ 15% in BaTi03, to be compared with 
- 5% in PbTi0 3 . 



B. Lattice instabilities 

As a first step towards understanding the unexpected 
behavior of PbTi03 at negative pressures, we investigate 
the zone-center lattice instabilities in a structure having 
the lattice constants of the optimized tetragonal struc- 
ture at the corresponding pressure, but without internal 
distortions (i.e., all atoms at centrosymmetric positions). 
For each structure, we compute the force-constant matrix 
at the r point using the frozen-phonon method. 

The tetragonal ferroelectric phase of PbTi03 belongs 
to the C\ v (PAmm) space group. At the T point, 
the vibrational representation is spanned by two one- 
dimensional irreducible representations A\ and B\ and 
one two-dimensional representation E, of which there 
are 4, 1, and 5 copies, respectively, so that the force- 
constant matrix is block diagonal. As we are interested 
in the modes producing polarization along z, we only 
have to diagonalize the 4x4 block of the Ai subspace. 
The pure translational mode is discarded, and the eigen- 
values (k) of the three remaining A\ modes are plotted 
versus pressure in Fig. 0] Negative values correspond to 
unstable modes. As expected, the ferroelectric soft mode 
is already unstable at zero pressure, and it becomes even 
more so as the crossover negative pressure region around 
p c = —4.8 GPa is crossed because of the cell-volume en- 
hancement that takes place there. However, we also find 
that a second mode becomes soft and that its frequency 
crosses through zero at a pressure that is close to the 
same p c . At first sight, this concurrence might be taken 
as a hint of some connection between the second mode 
crossing and the anomalous "kinky" behaviors. However, 
as will be explained in Sec. IV, we think that if there is 
such a connection, it is more likely that the rapid ex- 
pansion of the c lattice constant causes the second mode 
crossing, and not vice versa. 

The corresponding eigenvectors are also very sensitive 
to pressure changes. Table [I] summarizes the eigenval- 
ues k and the eigenvectors of the original ferroelectric 
soft mode and the second soft mode at three different 
external pressures: zero, p c , and a pressure higher than 
p c - As can be seen, the original soft mode shows the 
displacement pattern that is typical of ferroelectric per- 
ovskites, with the cations moving in opposition to the 
oxygen octahedra. As the negative pressure increases, 
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FIG. 4: Eigenvalues k of the force-constant matrix in the 
centrosymmetric tetragonal structure with lattice constants 
of the optimized tetragonal structure at the specified pressure. 
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FIG. 5: Comparison of lattice parameters c and a, in a.u., 
as a function of absolute value of stress obtained for tetrag- 
onal PbTiOs by applying negative biaxial stress or positive 
isotropic stress (negative pressure). 



the character of the mode evolves to one dominated by 
the long-short alternation of the bonds in the Ti-0 bonds 
along z. On the other hand, the pattern of the second soft 
mode consists mainly of the Pb-C>3 plane moving against 
the Ti-Oi plane. At zero pressure the opposing cation 
displacements are dominant; negative pressure leads to 
an increase in the Oi,2 displacements and a decrease in 
the Ti displacement. 



C. Applied biaxial stress 

As mentioned in Sec. Ill A, the effect of negative 
isotropic pressure on the cell shape is to stretch the unit 
cell along z and slightly squeeze it in-plane. Here we con- 
sider what happens if instead we apply an in-plane stress 
only (a xx =a yy and cr zz —0). This is effectively what hap- 
pens in the case of epitaxial growth on a substrate of 
slightly different lattice constant H In the calculation, we 
constrain the a lattice parameter and allow c and all the 
ionic positions in the [001] direction to relax while pre- 
serving the tetragonal space-group symmetry; the biaxial 



TABLE I: Comparison of eigenvalues (k) and eigenvectors 
of the original and second soft modes for different pressures. 
Eigenvectors (Pb, Ti, Oi=02, O3) are normalized to unity. 



p (GPa) k (a.u.) 



Eigenvector 



Original 



Second 



-0.187 
-4.56 -0.308 
-6.48 -0.383 





-4.56 
-6.48 



0.190 
0.056 
-0.116 



0.520 
0.217 
0.093 

0.713 
0.824 
0.814 



0.573 
0.752 
0.755 

-0.682 
-0.321 
-0.128 



-0.388 -0.317 

-0.211 -0.546 

-0.108 -0.631 

-0.076 0.121 

-0.317 0.130 

-0.394 0.102 



stress needed to support this structure is then computed 
directly within ABIIMIT. 

The results are shown in Fig. [SI where the lattice pa- 
rameters are plotted for expansive (positive) isotropic 
stress (negative pressure) and for compressive (negative) 
in-plane stress. This is the appropriate comparison since 
both have the effect of enlarging the c lattice constant and 
thus potentially triggering the abrupt structural change. 
Indeed, we find that the biaxial stress does generate the 
same kind of anomalous c-axis behavior. The enlarge- 
ment in c occurs at a somewhat larger magnitudes of 
the stress, and produces a somewhat more modest c-axis 
expansion, than for the isotropic-pressure case. 

It should be emphasized, however, that the abrupt 
variation only appears when c and a are plotted vs. ap- 
plied biaxial stress; simply plotting c vs. in-plane strain 
(i.e., vs. a) does not reveal anomalous behavior. More- 
over, in order to obtain the same large c values that result 
near the kink in the negative-pressure case (indicated by 
p c in the figure) , we find that the in-plane lattice param- 
eters would have to be compressed by ~ —3.2%, which 
is much larger than the typical compression that can be 
attained by growth of PbTiOs on typical substrates such 

as SrTiOs ( 1.4%)±£ or NdGa0 3 ( 1%). 19 Thus, 

it does not appear that the predicted behavior can be 
observed in epitaxially strained films. Perhaps it may be 
possible to find some other way to generate a sufficiently 
large biaxial compressive stress so that this effect can be 
observed. 

In any case, we think it is interesting that the applica- 
tion of compressive biaxial stress can generate the same 
kind of anomalous structural behavior as the application 
of negative isotropic pressure, at least for PbTiOs. This 
result tends to support the speculation that some other 
kinds of variation (e.g., chemical substitution) might also 
be capable of driving a similar anomalous behavior. 
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D. Applied electric field 

Another means of producing an effective tensile stress 
along z is the application of an electric field in the 
same direction. The induced polarization should cou- 
ple strongly to the tetragonal strain. Therefore, we 
made a preliminary investigation of whether a strong ap- 
plied electric field might be capable of inducing a similar 
anomalous enhancement of tetragonality. The structural 
parameters in applied electric fields can be estimated us- 
ing an approximate scheme in which the ab-initio force on 
each atom (computed in zero electric field) is augmented 
by the product of the dynamical effective charge tensor 
for that atom and the applied electric field vectori&211 
Based on our preliminary results, but more specifically 
on the results plotted in Fig. 8 of Ref. H, it appears 
that 7/3 and r\\ do not show any abrupt change up to 
6 x 10 3 kV/cm. Thus, it seems unlikely that the ap- 
plication of an electric field can cause the same kind of 
anomalous behavior as the application of negative pres- 
sure. A more careful study would be desirable, checking 
whether the above approximations might be oversimpli- 
fying the treatment of some nonlinearities, but this is left 
for a future investigation. 



IV. PHENOMENOLOGICAL MODEL 
DESCRIPTION 

In order to explore the origins of the anomalous behav- 
ior that has emerged from our calculations in the vicinity 
of the crossover negative pressure p c , we turn now to an 
attempt to model this behavior phenomenologically in 
terms of a reduced set of relevant degrees of freedom. 
We consider models similar to those that underlie the 
effective-Hamiltonian schemejSiSi in which the total en- 
ergy is Taylor-expanded, in soft-mode and strain vari- 
ables, about a reference cubic phase. The effect of the 
external hydrostatic pressure p is included by minimizing 
an enthalpy that includes a pV term in addition to the 
energy. 

The first step in the construction of the model is to 
determine an appropriate set of degrees of freedom to 
be included. For practical reasons, this set should be 
as small as possible, and the definition of the degrees 
of freedom should not depend explicitly on pressure. 
Starting from a reference cubic structure, the fully 
relaxed structure at a given pressure can be separated 
into two parts: a homogeneous strain (leaving the atoms 
undistorted from their high-symmetry positions) and an 
"internal strain" (i.e., internal atomic displacements). 
To describe the first part, we make the simplest possible 
choice. As the unit-cell volume V increases monotoni- 
cally with negative pressure, it specifies the strain state 
uniquely This choice has the added convenience that 
it is V that directly couples to the pressure (it should 
be noted that the cell at given V will in general be 
tetragonal and its value of c and a are thus functions 
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FIG. 6: Projection of the relaxed atomic displacements along 
z at the specified pressure onto the three zone-center 

modes of the zero-pressure cubic structure (\<Pu), \<fiv) and 
\(fiw)), given in units of c. 



of V). Next, we need to find an appropriate basis able 
to capture the important internal strains at all relevant 
pressures. We first consider the entire set of zone-center 
modes of the cubic structure at zero pressure, and 
then attempt to determine which of these are most 
relevant for spanning the observed configurations in the 
pressure range of interest. Since the structure remains 
tetragonal, the three non-trivial A\ modes form the 
starting point for this investigation. Let u, w and v 
denote the amplitudes of the softest, intermediate, and 
hardest mode of the force-constant matrix, with eigen- 
vectors |<p u ) = |0.66,0.39,-0.43,-0.43,-0.21), 
\<p w ) = |0.60,-0.79,0.10,0.10,-0.01) and 

\<p v ) = |0.06,0.15,0.33,0.33,-0.87), respectively. 
Then for each pressure, the optimized displacements 
shown in Fig. ^c) are projected onto these three modes 
to obtain the amplitudes plotted in Fig. As can 
be seen from the figure, the largest contribution at all 
pressures is clearly given by the soft mode u, while 
the other two appear less important. The contribution 
of the second mode is almost independent of pressure 
and can therefore be neglected. On the contrary, the 
highest-mode contribution is quite sensitive to pressure, 
increasing notably at p c . The need to include at least 
one additional mode in the model subspace is clear 
from Figure 1(c), as a single mode cannot capture the 
change in character above p c . The figure suggests that 
the highest-mode contribution v is next in importance 
after the soft mode u; the fact that including the v 
mode yields a qualitatively correct description of the 
distortions at all relevant pressures will be confirmed 
below. 

To develop the model, we consider the configuration 
subspace defined by three independent variables: the 
unit-cell volume V, the soft-mode amplitude u, and 
the hardest-mode amplitude v. At each volume V, we 
start by considering only the mode u (setting v=0), and 
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TABLE II: Parameters of the model total-energy expansion. 
All parameters are in atomic units. 



V 


A 


4.4668 


Ai 


-0.0124 




A 3 


2.37xl0" 7 


Ai 


-2.17xl0" 10 


u 




0.9264 




-0.0042 




Cau 


0.7095 




-0.0030 


V 


Bi v 


0.6700 


B[ v 


-1.64xl0~ a 




B 2 v 


-0.0851 


B' 2v 


5.16xl0" 4 



A 2 



-5.69x10" 



r" 
r" 



4.65x10 
3.18x10 



Taylor-expand in u as 

E(V, u) = E + c (V) + c 2u (V) u 2 + c iu (V) u 4 , (1) 

where Eq is the total energy at V = Vq and u = v = 0, 
so that co(Vb) = 0, and terms of order u 6 and higher are 
dropped. The coefficients cq, c 2u and C4„ are obtained 
by fitting to ab-initio total energies computed for con- 
figurations corresponding to various values of m, and the 
equilibrium value for v = and fixed V is obtained by 
minimizing Eq. l[T|l. leading to 



< q {V) = -- 



1 C 2u (V) 



2 c 4m (V) 



and 



E(V, u eq (V)) = E + c (V) 



4 c 4m (V) 



(2) 



(3) 



Then, starting at the configuration corresponding to 
(u, v) = (u eq (V),0), we compute ab-initio total energies 
for configurations corresponding to various values of v, 
holding V and u — u eq (V) fixed, and fit the result to 
quadratic order in v as 

E(V, u eq (V), v) = E(V, u eq {V)) + b lv (V) v + b 2v (V) v 2 . 

(4) 

Minimizing the total energy with respect to v (at fixed 
V and u = u eq ), the equilibrium amplitude of this mode 
is obtained as 



v eq (V) = -- 



1 b lv (V) 



and the total energy is 
E = E Q + c (V) 



2 b 2v (V) 



l c 2 u (V) l b 2 lv (V) 
4 c 4u {V) 4 b 2v [V) 



(5) 



(6) 



To facilitate the modeling of the dependence of these 
results on volume V, we fit the functions cq, c 2u , C4„, b\ v , 
and b 2v as polynomials in V: 



Co (V) = A Q + A\V + A 2 V 2 



A,V 3 



A A V A 



c 2u {v) = c 2u + c 2u v + c'i u v 2 



cau (v) = c iu + c' 4u v + c'ly 2 



(7) 



(8) 



(9) 
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FIG. 7: Pressure computed at different volumes by using the 
phenomenological model (full line) compared with the first- 
principles results (circles). 



bi v (V) - B lv + B' lv V , 



b 2v (V) = B 2v 



B' 2v V 



(10) 



(11) 



The expansion parameters were determined by fitting to 
the results of total-energy calculations in the interval of 
volumes from 400 to 480 a.u. 3 , near V = 393.99 a.u. 3 . 
The resulting parameters of the model are reported in 
Table HH 

We then explored the behavior of this model by deter- 
mining the equilibrium structural parameters as a func- 
tion of pressure p, minimizing the enthalpy H = E +pV. 
The result is shown by the solid line in Fig. We 
find that there is a pressure interval from about —4.6 
to —6 GPa in which two local minima compete, with 
a first-order isostructural transition between minima at 
p c ~ —5.4 GPa. The locations of the secondary min- 
imum and the saddle point are indicated by the dotted 
line in Fig.0 obtained by evaluating the thermodynamic 
relation p= -dE(V)/dV function of V. 

Comparing the behavior of the model with the first- 
principles results, shown as solid circles in Fig. \7\ we see 
that the model correctly reproduces the existence of an 
abrupt variation of structural parameters with negative 
pressure. In fact, the model even goes too far, exhibiting 
a true first-order transition where there is none in the 
first-principles results. The model also slightly overesti- 
mates, by about 10-15%, the magnitude of the negative 
pressure at which the abrupt change occurs, and under- 
estimates the pressure dependence of the volume at the 
highest negative pressures. 

Even with these minor discrepancies, the model is very 
useful in generating a clearer picture of the observed be- 
havior. In particular, the fact that a first-order transition 
occurs in the model supports the idea that the anomalous 
enhancement of tetragonality in the first-principles calcu- 
lations results from proximity to a phase transition. To 
make this idea more precise, we consider small variations 
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FIG. 8: Calculated P — V phase diagram for different C'i u 
parameters (upper panel). The original model is marked as 
case (1). C'l u parameter-critical pressure (p c ) diagram (bot- 
tom panel). For values bigger than C'l u at the critical point 
there is no phase transition. 



of one of the model parameters, C'l u , keeping all other 
parameters fixed, and study the transition behavior as 
a function of pressure. The upper panel of Fig. [S] shows 
the P-V diagram for four different values oiC'l u , starting 
from the original model marked as case (1) and increas- 
ing C'l u incrementally to case (4). As can be seen, for the 
smaller values of C'( u there is a first-order phase transi- 
tion from a low-c/a to a high-c/et phase at a transition 
pressure p c (cases (1) and (2)). The volume discontinu- 
ity decreases with increasing C" u and vanishes at a triple 
point (case (3)). In other words, a critical value of the 
C'l u parameter exists above which there is no distinction 
between two phases (see, e.g., case (4)). The relation be- 
tween C'l u and the transition pressure is presented in the 
bottom panel of Figure [H] the phase boundary defined 
in this way terminates at the triple point identified by 
pressure p* . 

Thus, Fig. [S] demonstrates that with increasing nega- 
tive pressure, PbTiC>3 passes very close to a triple point. 
In fact, it passes so close that the errors introduced by the 
truncations and simplifications of our model, described 
earlier in this section, are sufficient to cause the model 
system to exhibit a true first-order phase transition. In- 
deed, to shift the model system into the transition-free 
region of parameter space, it suffices to change C'l u by 
only ~ 0.15%. A more accurate model containing ad- 



FIG. 9: Configurational space described by u eq (upper panel) 
and by u eq and v eq (bottom panel) compared with the 
minimum-energy structures (dotted lines). 



ditional fitting parameters (e.g., an independent treat- 
ment of lattice constants c and a, inclusion of mode w, 
and/or a more careful treatment of cross terms between 
modes u and v) would presumably be sufficient to repro- 
duce the observed "kinky" behavior without the spurious 
first-order transition. 

We now turn to an examination of the role of the addi- 
tional mode v in the anomalous negative-pressure behav- 
ior of PbTiC>3. The effects of holding v—0 can be easily 
obtained by minimizing the enthalpy E+pV with E from 
Eq. ©. The V vs. p curve still exhibits the anomalous 
negative-pressure behavior, although the transition re- 
gion is considerably shifted (to around —12 GPa) and 
there is no first-order transition in this case (similar to 
the first-principles result). Thus, inclusion of the ad- 
ditional mode v is not absolutely necessary to produce 
the anomalous negative-pressure behavior. On the other 
hand, neglecting this mode leads to a substantial quan- 
titative error in the transition pressure. Moreover, as we 
have already seen in Fig. [He) an d Fig. [S] the inclusion 
of the mode v is needed to correctly capture the vari- 
ation of the atomic displacements with pressure. This 
point is emphasized in Fig. In the upper panel, the 
ratios of the displacements are fixed with an overall am- 
plitude proportional to u (shown with symbols), so that 
it is impossible to reproduce the crossing in the oxygen 
displacement pattern around p c . In the lower panel, the 
inclusion of the v mode yields a qualitative improvement, 
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with a crossing close to p c . The discrepancy between the 
model and the first-principles results (shown with dotted 
lines) grows much more noticeable above p c . This most 
likely reflects the contribution of the neglected modes 
and couplings to the strong pressure dependence above 
Pc which is underestimated by the model. 

Finally, we analyze the relation between the modes of 
the centrosymmetric tetragonal structure (dependent on 
pressure) and the modes u, v and w, used in the construc- 
tion of our model. In our discussion of Fig.0]in Sec. Ill B, 
we noted that the tetragonal mode of intermediate hard- 
ness crosses from positive to negative K (i.e., from stable 
to unstable behavior) at a pressure very close to that at 
which the anomalous structural response occurs. Pro- 
jecting the pressure-dependent tetragonal modes on the 
cubic modes, we find that the tetragonal soft mode is 
mainly described by the mode u at all pressures, with a 
significant contribution of the mode w above p c , whereas 
the tetragonal mode of intermediate hardness is mainly 
described at lower pressures by the mode w, and then 
above p c by the mode v. Since the kinky behavior oc- 
curs even when both cubic-structure modes v and w are 
frozen to zero amplitude, we think that it is unlikely that 
the zero crossing of the second tetragonal-structure mode 
has a causal role in the anomalous structural behavior. 
Instead, it seems more likely that it simply results from 
the abrupt expansion of the c-axis lattice parameter. 

Thus, we believe that the model developed in this sec- 
tion, based on an expansion of the total energy with 
respect to the volume and two zone-center mode am- 
plitudes, is a reasonable compromise between simplicity 
and accuracy. It gives a good qualitative and even semi- 
quantitative description of the abrupt variation of the 
structural parameters near p c . The fact that it predicts 
a true first-order transition at p c , instead of a smooth 
but "kinky" cross-over, should be taken more as a sign 
of the proximity of the real system to a phase transi- 
tion than as a failure of the model. We are hopeful that 
this model will be useful for investigating and describing 
similar behavior in other perovskite systems. 

V. DISCUSSION AND SUMMARY 

We have shown that the application of negative pres- 
sure induces a large enhancement of tetragonal strain in 
tetragonal PbTiC>3. Specifically, in a window of pres- 
sure centered around a particular value p c , c/a abruptly 



increases, and all structural parameters exhibit a cor- 
responding "jump" that suggests proximity to a phase 
transition. 

To describe this unexpected behavior in PbTi03, 
we have generated a phenomenological model from our 
ab initio results, based on an expansion of the total en- 
ergy with respect to a carefully selected subset of the 
structural degrees of freedom. This has led to the identi- 
fication of a phase boundary, analogous to the liquid-gas 
phase boundary, in the model parameter space. With the 
parameters obtained from fitting to first-principles total 
energies of PbTiC>3, the model predicts not just a "kinky" 
behavior, but a true first-order transition at a critical 
pressure p c close to the pressure at which the anomaly 
occurs in the first-principles calculations. We found that 
a tiny variation of one parameter is sufficient to drive the 
model through a triple point and bring it to a transition- 
free portion of parameter space where the "kinky" be- 
havior is qualitatively well-reproduced. It might become 
possible to realize such variations of effective model pa- 
rameters by some combination of "external fields" such 
as chemical substitution, temperature, epitaxial stress, 
or homogeneous electric field. If so, one might develop 
experimental techniques for turning the transition off or 
on, tuning the associated transition properties, and per- 
haps even exploring the region of the triple point itself. 

Finally, we speculate that this "anomalous" behavior 
of the structural parameters of tetragonal PbTiC>3 could 
be a more general feature in perovskite oxides. For ex- 
ample, we have found that the same anomalous behav- 
ior is also present in BaTiC>3 with imposed tetragonal 
symmetry (though at much higher negative pressure), so 
that a similar phase boundary could be explored in this 
portion of perovskite parameter space. The possibility of 
analogous transitions in other ABO3 pcrovskites is under 
investigation. This could present tantalizing opportuni- 
ties for designing perovskite-based materials with large 
and controllable strain variations. 



Acknowledgments 

We would like to thank Na Sai for providing the ini- 
tial suggestion for the direction of this work. This work 
was supported by the Center for Piezoelectrics by Design 
(CPD) under ONR Grant N00014-01-1-0365. Computa- 
tional facilities for the work were also provided by the 
CPD. We thank Morrel Cohen for useful discussions. 



1 S.-E. Park, T. R. Shrout, J. Appl. Phys. 82 1804 (1997). S.- 
E. Park, W. Hackenberger, Curr. Opin. Solid State Mater 
Sci. 6, 11 (2002). 

2 B. Noheda, D. E. Cox , G. Shirane, S.-E. Park, L. E. Cross, 
Z. Zhong, Phys. Rev. Lett. 86 3891 (2001). 

3 B. Noheda, D. E. Cox , G. Shirane, J. Gao, Z. G. Ye, Phys. 
Rev. B 66 054104 (2002). 



4 R. E. Cohen, Nature 258, 136 (1992). 

5 A. Garcia and D. Vanderbilt, Phys. Rev. B 54, 3817 

(1996) . 

6 U. V. Waghmare and K. M. Rabe, Phys. Rev. B 55, 6161 

(1997) . 

7 G. Saghi-Szabo, R. E. Cohen and H. Krakauer, Phys. Rev. 
Lett. 80, 4321 (1998). 



9 



Na Sai, Karin M. Rabe and David Vanderbilt, Phys. Rev. 
B 66, 104 108 (2002). 
9 ABINIT is a common project of the Universite Catholique 
de Louvain, Corning Incorporated, and other contribu- 
tors (http://www.abinit.org). X. Gonze, J.-M. Beuken, 
R. Caracas, F. Detraux, M. Fuchs, G.-M. Rignanese, L. 
Sindic, M. Verstraete, G. Zerah, F. Jollet, M. Torrent, A. 
Roy, M. Mikami, Ph. Ghosez, J.-Y. Raty, D.C. Allan, Corn- 
put. Mater. Sci. 25, 478-492 (2002). 

10 J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13 244 (1992). 

11 D. Ceperley, Phys. Rev. B 18, 3126 (1978); D. M. Ceperley 
and B. J. Alder Phys. Rev. Lett. 45, 566 (1980). 

12 M. Teter, Phys. Rev. B 48, 5031 (1993). 

13 G. Shirane, R. Pepinsky and B. C. Frozer, Acta Crystal- 
logr. 9, 131, (1956). 

14 R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 



1651 (1993). 

15 W. Zhong, R. D. King-Smith and D. Vanderbilt, Phys. 
Rev. Lett. 72, 3618 (1994). 

16 Ph. Ghosez, J. P. Michenaud, X. Gonze, Phys. Rev. B 58, 
6224 (1998). 

17 N. A. Pertsev, A. G. Zembilgotov, A. K. Tagantsev, Phys. 
Rev. Lett. 80, 1988 (1998). 

18 Y.-F. Chen, L. Sun, T. Yu, J.-Y. Chen, N.-B. Ming, D.-S. 
Ding, L.-W. Wang, Appl. Phys. Lett. 67, 3503 (1995). 

19 L. Sun, Y.-F. Chen, L. He, Ch.-Z. Ge, D.-S. Ding, T. Yu, 
M.-S. Zhang, N.-B. Ming, Phys. Rev. B 55, 12 218 (1997). 

20 H. Fu and L. Bellaiche, private communication. 

21 W. Zhong, D. Vanderbilt and K. M. Rabe, Phys. Rev. Lett. 
73, 1861 (1994). 



